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ABSTRACT 


w 

’  A  1-d  electrostatic  and  relativistic  particle  code  has  been 
built  in  order  to  simulate  in  a  self-consistent  way  the  time-dependent 
behavior  of  collective  ion  acceleration  produced  by  linear  electron 
beams  injected  into  evacuated  drift  tubes.  The  simulation  results 
agree  with  existing  experimental  and  theoretical  results  by  showing 
good  current  propagation  and  high  ion  energy  gain  in  the  presence  of 
a  well-localized  plasma  source  at  the  e-beam  injection  boundary.  More¬ 
over,  the  simulations  display  an  unexpected  phenomenon,  the  appearance 

S 

of  short-lived  regions  of  positive  electrostatic  potential  ('Virtual 

e 

anodes")  which  together  with  the  usual  virtual  cathodes,  contribute 
to  the  ion  acceleration  mechanism.  In  the  absence  of  a  suitable 
analytical  theory,  it  seems  to  us  that  the  present  computer  model 
is  an  indispensable  tool  for  reliable  explanation  and  optimization 
of  present  and  future  collective  ion  acceleration  experiments.  Use 
of  our  model  in  laser-plasma  experiments  is  suggested.  ,\ 


I.  INTRODUCTION 


In  recent  years,  heavy  ion  beams  have  been  considered  as  possible  drivers 
in  an  inertial  confinement  fusion  scheme,  because  of  their  good  deposition  char¬ 
acteristics  and  because  of  the  already  existing  advanced  accelerator  technology, 
capable  of  high  repetition  rates  [1].  However,  a  serious  drawback  is  the  high 
cost  of  these  accelerators.  A  promising  and  possibly  less  expensive  approach 
is  the  use  of  collective  ion  acceleration  (CIA)  as  a  source  for  heavy  ion  beams. 
In  a  typical  CIA  experiment  a  relativistic  electron  beam  is  injected  into  a  drift 
tube  in  which  an  ion  source  is  present.  The  ions  which  initially  have  negligible 
kinetic  energy,  are  eventually  accelerated  by  the  collective  electrostatic  fields 
generated  by  the  incoming  electron  beam.  Basically,  in  this  process  long  e-beam 
pulses  are  transformed  into  short  energetic  ion  pulses. 

The  collective  acceleration  of  protons  and  heavier  ions  by  linear  electron 
beams  injected  into  initially  evacuated  drift  tubes,  has  been  investigated  at 
several  laboratories  [2-5].  The  ion  source  is  produced  either  from  an  insulating 
anode  or  from  a  puffed-in  gas,  by  an  e-beam  prepulse  or  possibly  by  laser  irra¬ 
diation.  Reported  experiments  indicate  that  effective  e-beam  propagation  and 
effective  collective  acceleration  occur  when  the  plasma  source  is  initially  well- 
confined  to  the  anode.  Collectively  accelerated  protons  have  been  observed  in 
these  experiments  with  energies  of  up  to  20  times  the  beam  electron  energy,  and 
heavier  ions  (C,  N,  0,  F,  Cl)  have  been  accelerated  to  energies  of  several  MeV 
per  nucleon. 

In  order  to  explain  the  CIA  mechanism,  it  has  been  recently  suggested  |.\M 
that  some  ions  from  the  dense,  anode-localized  plasma  source  are  trapped  and 
accelerated  by  the  electrostatic  fields  associated  with  a  moving  potential  well 
of  a  depth  comparable  to  the  injected  e-beam  energy.  An  explicit  relationship 
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is  given  for  the  fastest  ion  energy  as  function  of  time,  injected  energy  and 
current,  ion  mass  and  charge  state,  and  drift  tube  length.  The  effect  of  the 
ion  charge  on  the  potential  distribution  is  not  taken  into  account.  However, 
experiments  indicate  a  strong  coupling  between  electrons  and  ions  during  their 
propagation  towards  the  downstream  end  of  the  drift  tube  (coherent  movement 
is  observed) .  This  means  that  both  electrons  and  ions  may  have  comparable 
effects  on  the  development  of  the  electrostatic  potential  well. 

Therefore,  in  spite  of  the  success  of  the  above  theory  in  giving  some  crude 
evaluation  of  the  maximum  achievable  ion  energies,  there  is  a  critical  need  for 
a  time-dependent  model,  in  which  ions  and  electrons  are  included  in  a  self-con¬ 
sistent  manner.  In  order  to  investigate  the  properties  of  the  fully  nonlinear 
ion  acceleration  mechanism  we  have  built  a  one-dimensional  relativistic 
particle  simulation  code,  which  includes  the  dynamic  equations  for  the 
electrons  and  ions  simultaneously.  This  should  allow  a  more  detailed  assess¬ 
ment  of  the  ion  energy  gain.  Besides  being  the  only  means  for  a  self-consistent 
time-dependent  study  of  CIA,  computer  simulation  allows  for  a  multitude  of  diag¬ 
nostics  and  a  strict  control  of  the  system  parameters.  Once  the  acceleration 
mechanism  is  understood,  we  can  determine  the  optimum  system  parameters  for 
experiments  planned  at  several  laboratories  [2-4],  in  order  to  accelerate  ions 
to  maximum  energy  and  current. 

In  Section  II,  we  describe  the  technical  details  of  our  simulation  code. 

In  Sec.  Ill,  we  present  preliminary  simulation  results  and  compare  them  with 
theory  and  experiments.  In  Sec.  IV,  we  outline  the  main  conclusions  of  this 
study. 


II.  THE  SIMULATION  CODE 


A  sketch  of  the  simulated  physical  system  is  presented  in  Fig.  1.  A  relativ¬ 
istic  electron  beam,  usually  produced  by  a  diode  is  injected  into  an  initially 
evacuated  drift  tube,  forming  a  virtual  cathode  (defined  as  the  minimum  of  the 
potential  well).  If  an  unlimited,  infinitely  thin  plasma  layer  is  available  at 
the  injection  site,  ions  will  be  emitted  from  the  plasma  source  in  a  space-charge 
limited  manner  and  will  be  accelerated  by  the  virtual  cathode  towards  the  down¬ 
stream  end  of  the  drift  tube.  Whenever  the  electric  field  at  the  injection  boundary 
becomes  negative,  electrons  will  be  emitted  from  the  plasma  source,  instead  of 
ions.  The  electrons  and  the  ions  in  the  drift  tube  are  free  to  move  under  the 
effect  of  the  self-consistent  collective  electric  field  generated  by  the  net  charge 
density  distribution.  The  potential  difference  across  the  tube  is  assumed  to  be 
always  zero,  and  the  charge  emission  occurs  due  to  the  presence  of  the  potential 
well  generated  by  the  incoming  electron  beam.  The  injected  electrons  have  all 
the  same  velocity,  corresponding  to  their  relativistic  energy.  The  emitted  ions 
or  electrons  start  with  negligible  kinetic  energies.  In  the  present  code  we 
assume  the  injected  current  and  energy  to  be  constant,  although  any  time-dependence 
can  be  accomodated  in  our  model.  Both  electrons  and  ions,  emitted  or  injected 
into  the  drift  tube  are  absorbed  by  the  downstream  side  of  the  tube,  once  they 
reach  it.  This  side  of  the  drift  tube  is  not  assumed  to  emit  particles,  although 
this  may  happen  in  reality.  Therefore,  the  electric  field  is  strictly  zero  only 
at  the  injection  boundary  and  at  the  minima  of  the  potential  well.  Usually,  the 
source  of  the  injected  electron  beam  is  a  relativistic  diode,  whose  anode  is  i- 
dentical  with  the  injection  boundary  of  the  drift  tube.  Our  code  simulates  the 
drift  tube  part  of  the  system,  with  the  e-beam  current  and  energy  externally  given. 

The  physical  parameters  in  the  simulation  are  as  follows:  the  e-beam  energy. 


VQ(t);  the  e-beam  current  density,  Jg(t);  the  drift  tube  length,  d;  and  the  charge 

to  mass  ratio  of  the  ions,  q^/m^.  The  simulation  parameters  are:  the  space  step, 

A z;  the  time  step.  At;  the  number  of  particle  records  (half  electrons,  half  ions), 

NREC;  the  number  of  particles  per  record,  NBATCH;  the  number  of  cells  in  the  system, 

NZ  (d  =  NZ  •  Az);  the  number  of  particles  to  be  emitted  at  each  time  step,  NEMIT. 
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Our  PIC  code  is  1  y  dimensional  (z,  vz,  vr),  electrostatic  and  relativistic.  A 

flow  chart  of  the  simulation  code  is  presented  in  Fig.  2.  Initially,  the  code  starts 

2 

with  a  vacuum  drift  tube:  p  (z)  =  0.  The  code  solves  the  poisson  equation  V  <f>(z)  = 
-4irp(z),  where  <p (z)  is  the  potential  distribution,  and  p  (z)  is  the  net  charge  dis¬ 
tribution  in  the  system.  The  boundary  conditions  are:  cp (Al)  =  tf>(A2)  =  0,  where 
A1  and  A2  are  the  injection  and  the  downstream  boundaries,  respectively.  After 
finding  <f>(z),  the  electric  field  Ez(z)  is  calculated  from: 

Ez(z)  =  U(z)  -  <j>(z  +  Az)]/Az  . 

This  gives  Ez(Al),  the  electric  field  at  the  injection  boundary,  before  emission. 

Then  we  use  the  gaussian  law  for  space  charge  limited  emission  at  Al,  in  order  to 

make  the  electric  field  equal  to  zero,  at  the  emission  boundary.  If  E  (Al)  is 

z 

positive,  ions  are  emitted  at  Al,  according  to  the  relation: 

Qjnit(Al)  =  Ez(A1)/4ttA  , 

where  is  the  surface  charge  density  of  the  emitted  ions,  and  A(  =  1  -  0.5/NZ) 

is  a  correction  factor  aimed  at  matching  the  boundary  conditions  to  the  finite 
grid  system.  If  Ez(Al)  is  negative,  electrons  will  be  emitted  at  Al,  (in  addition 
to  the  injected  electrons),  given  by  the  same  expression: 


Q  i  (Al)  =  E  (Al)  / 4ttA 


The  corresponding  current  densities  are  defined  as: 

J^"lt(Al)  =  Q^Tnit(Al)/At 


jJJlt:(Al)  =  QjJlt:(Al)/At 

sbs  dbs 

The  other  current  density  components  in  the  system  are:  (Al) ,  J  (A2), 

aKo  aKo 

(Al),  (A2),  where  el  =  electrons,  i  =  ions,  abs  =  absorbed.  No  currents 

are  assumed  to  be  emitted  at  A2.  Q^-1  is  the  surface  charge  density  of  the  in¬ 

jected  electrons  defined  as: 

Qelj  =  ~JenjAt/cMKS, 

where  CMKS  converts  amps  to  CGS  units  if  J*n;I  is  given  in  A/cm2;  is  the  injected 

electron  current  density. 

Next,  the  electrostatic  potential  is  modified  to  take  into  account  the  effect 
of  the  emitted  and  injected  charges  at  Al.  This  is  done  by  using  a  Green-like 
function  which  represents  the  correction  to  the  potential  due  to  one  unit  of  sur¬ 
face  charge  density  at  Al.  This  function,  A<j> e ,  does  not  alter  the  potential 
boundary  values.  It  only  makes  the  electric  field  equal  to  zero  at  the  emitting 
electrode,  Al.  This  correcting  function  is  given  by  [8]: 


A  e  (  z  ) 


=  -2ttAz  (z-d) /d ,  +Az/2  <  z  <  d 


=  -4ir(Az/2-d)/d,  -Az/2  <  z  <  Az/2 


The  corrected  electrostatic  potential  function  will  then  be: 


tC(z)  =  *(z)  +  [Qjmlt (Al)  +  Qe^nit(Al)  +  Q*"j(Al)]  •  A<|>e(z) 


The  corrected  electric  field  Ez(z)  is  therefore  given  by: 


E^(z)  =  (1/Az)  •  [$C(z)  -  <f>C(z  +  Az)], 
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which  satisfies:  EC(A1)  =  0. 

z 

The  particles  are  then  actually  emitted  into  the  system.  The  injected 
electrons  have  all  the  same  longitudinal  velocity,  (no  transverse  momentum), 
given  by: 

v*  ■  c-/2  V0<1  +  0.5  »„)/(!  +  Vo)> 


where  V  =  V  /m  c 
o  o  e 


and  initially  they  are  placed  according  to: 


z  =  0.5  •  random  (0,1)  +  At  •  v6/2Az 

z 

p 

The  electric  field  Ez(z)  is  then  used  to  push  the  particles  (relativistically 
for  electrons  or  positrons,  nonrelativistically  for  ions)  to  their  new  positions, 
during  the  time  step  At.  The  new  charge  distribution  p(z)  is  calculated  by  the 
usual  grid-weighting  method,  and  finally  the  code  proceeds  to  the  next  time  step. 

Many  diagnostics  are  built  in  the  code  for  analyzing  the  electrostatic 
field  and  potential,  the  various  current  components  and  the  energy  spectra 
of  the  particles.  Various  options  are  available  within  this  code,  as  follows: 

(a)  Ions  can  be  emitted  in  various  modes,  corresponding  to  the  experi¬ 
mental  designs.  For  example,  besides  the  space-charge  limited  emission  men¬ 
tioned  above,  ions  can  be  emitted  from  a  finite  plasma  source  near  the  injection 
boundary,  until  the  source  is  depleted. 

(b)  A  delay  can  be  allowed  between  the  times  of  plasma  source  formation 
and  electron  beam  injection. 

(c)  The  potentials  at  A1  and  A2  may  not  be  equal.  In  fact,  experiments 

have  been  done  with  <J>(A1)  =  0  and  d>(A2)  =  -V  . 

o 

(d)  The  injected  electrons  need  not  be  mono-energetic,  but  can  have 

a  given  velocity  distribution  function  f(v  )  (e.g.  Maxwellian).  This  enables 
one  to  apply  the  code  to  the  analysis  of  fast  ions  produced  by  the  high  energy 
tail  of  the  electron  distribution  function  generated  in  laser  produced  plasmas. 
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The  choice  of  the  simulation  parameters  should  fulfill  some  basic  requirements. 
First,  the  time  step  should  be  small  enough  to  allow  rapid  time  variations  of  the 
physical  quantities.  For  example,  we  require  Avg  <  v™aXftc)  at  anytime,  where  vg 
is  the  electron  speed  and  c  is  the  light  speed.  Taking  the  maximum  possible  value 
for  the  electric  field  (when  the  voltage  between  A1  and  the  virtual  cathode  spans 
across  a  few  simulation  cells)  we  get:  At  <  (Az/Vq) (m ^c/e).  This  is  formally  similar 
to  a  Courant-type  condition,  Az  >  cAt,  in  which  the  fastest  electrons  are  not  allowed 
to  travel  over  more  than  one  cell  during  one  time-step  in  order  to  correctly  sample 
the  space-varying  field.  In  cases  with  high  injected  current  densities,  situations 
with  high  and  rapidly  space  and  time  varying  electric  fields  often  occur.  Thus, 
quite  small  time  steps  must  be  used  in  these  cases.  On  the  other  hand,  very  short  time 
steps  cause  prohibitively  long  computational  times  or  large  accumulation  errors. 

A  separate  condition  exists  on  Az,  to  ensure  that  the  system  can  accomodate 
the  highest  physically  attainable  currents  in  the  system.  This  condition  poses 
an  upper  limit  on  Az,  which  for  a  given  d,  determines  the  minimum  number  of  simu¬ 
lation  cells,  NZ,  to  be  used  (d  =  NZ  •  Az)  for  given  and  d. 

A  different  problem  exists  concerning  the  physical  parameters  of  the  simulations. 
In  all  our  simulations,  we  use  realistic  values  for  Jln  ,  V  and  d.  However,  for 
mi^me’  We  ^ave  chosen  a  value  of  25,  much  lower  than  in  the  experiments.  We  do 
this  in  order  to  reduce  the  computational  times  to  acceptable  levels.  A  study 
of  the  scaling-up  of  the  simulation  results  with  nu/me  and  other  parameters  is 
necessary.  From  preliminary  results  it  seems  that  only  the  time  scales  are 
affected,  while  the  final  saturation  values  of  current  components  and  ion  energies 
are  relatively  independent  of  itu/m^. 


A 


^ >  ... .....  -- 
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III.  SIMULATION  RESULTS 


As  a  first  test  case,  we  performed  a  simulation  with  only  the  electrons 

present.  The  simulation  parameters  were  the  following:  the  injected  e-beam 

2 

current  density,  Jg  =  10  KA/cm  ;  the  energy,  V  =  1  MeV;  At  =  1.5  psec;  Az  =  0.07  cm; 
NZ  =  128  (d  =  Az  •  NZ  =  9.1  cm).  Initially,  a  deep  potential  well  of  2-3  times 
Vq  is  formed  which  lasts  for  a  few  plasma  periods.  No  electron  current  propagates 
during  this  period.  Then,  the  potential  stabilizes  and  the  virtual  cathode  ampli¬ 
tude  (defined  as  the  minimum  of  the  potential  well)  oscillates  up  and  down  in  the 
steady  state  around  a  nearly  constant  position  of  %1  cm  (Figs.  3,4).  These  oscil¬ 
lations  have  been  suggested  as  possible  sources  for  microwaves.  Following  the 
stabilization  time,  there  is  almost  complete  reflection  of  the  electron  current. 

(no  more  than  IX  propagates  to  the  downstream  end  of  the  drift  tube).  The 
stabilization  time  of  the  potential  corresponds  roughly  to  one  crossing  time 
for  the  electron  in  the  drift  tube.  Our  results  are  in  complete  agreement  with 
both  theory  and  experiments  [2]  which  have  shown  that  when  no  ion  source 
is  present,  the  propagated  crrent  should  be  and  in  fact  is  always  below  the 
limiting  value. 

In  a  different  run,  with  the  same  parameters,  we  provide  an  unlimited  proton 
source  at  the  injection  boundary,  having  a  mass  ratio  of  m^/ mg  =  25.  The  emission 
is  assumed  to  be  space-charge  limited.  The  total  simulation  time  is  6  nsec,  more 
than  sufficient  for  the  system  to  achieve  a  dynamic  steady  state.  We  have  obtained 
both  expected  and  unexpected  results.  In  agreement  with  theory  and  experiments, 
most  of  the  electron  current  (^80%)  propagates  to  the  downstream  end  of  the  tube 
(Fig.  5)  and  ions  having  energies  several  times  the  injected  energy  V  ,  are  ob¬ 
served  (Fig.  6).  Also,  synchronous  propagation  of  the  electrons  and  the  fast  ions 
is  evident,  suggesting  a  strong  electron-ion  coupling. 


Contrary  to  expectations,  the  virtual  cathode  does  not  display  a  clear  motion 
downstream  the  drift  tube.  Rather,  it  has  an  irregular  shape,  changing  continuously 


with  time.  Obviously,  the  most  interesting  results  are  those  which  have  not 
been  predicted  by  any  existing  analytical  theory.  First,  the  whole  ion  energy 
spectrum  (the  number  of  ions  in  different  energy  ranges),  has  been  obtained,  as 
a  function  of  time  (Fig.  6).  Second,  an  unexpected  feature,  the  formation  of 
"virtual  anodes"  is  observed,  identified  as  short-lived  regions  of  positive 
electrostatic  potential  (Fig.  7).  They  are  due  to  the  deceleration  and  bunching 
of  ions  in  the  region  between  the  virtual  cathode  and  the  downstream  end  of  the 
tube.  A  virtual  anode  is  capable  of  further  accelerating  the  ions  towards  the 
downstream  end  of  the  drift  tube.  It  also  accelerates  the  incoming  electrons. 

This  suggests  the  following  mechanism  for  collective  ion  acceleration:  first  the 
incoming  electrons  are  slowed  down  and  transfer  their  energy  to  the  initially 
slow  ions,  in  the  energy  range  0-V  .  Then,  those  ions  form  a  "virtual  anode",  which 
accelerates  the  electrons.  The  virtual  anode  and  the  accelerated  electrons  further 
accelerate  a  smaller  number  of  ions.  This  process  can  repeat  itself,  given  a  long 
enough  e-beam  pulse  duration  and  a  long  enough  drift  tube.  This  picture  is  further 
strengthened  by  the  ion  energy  spectra  and  their  time  derivatives  (Figs.  6 
and  8).  First,  one  can  see  that  the  number  of  slow  ions  (0-Vo)  is  continuously 
increased  by  the  incoming  e-beam.  Then  the  slower  ions  gradually  transfer 
their  energy  to  faster  ions.  More  details  can  be  found  in  Ref.  7.  Obviously 
with  these  and  other  diagnostics,  our  particle  simulation  code  is  capable  of 
following  in  detail  the  CIA  mechanisms. 


IV.  CONCLUSIONS 


Collective  ion  acceleration  (CIA)  is  a  complex  phenomenon.  We  have  built 
a  computer  simulation  model  including  the  dynamics  of  electrons  and  ions  simul¬ 
taneously  which  can  be  used  to  optimize  the  ion  energy  gain  in  present  CIA  experi¬ 
ments.  A  thorough  parametric  investigation  of  CIA  is  underway  as  a  function  of 
v  ,  d  and  other  physical  quantities.  With  our  code,  we  have  been  able  to 
identify  an  unexpected  phenomena,  the  formation  of  short-lived  "virtual  anodes", 
capable  of  accelerating  ions  and  electrons  to  high  energies.  Preliminary  simu¬ 
lations  show  good  qualitative  agreement  with  experimental  results. 

In  the  absence  of  a  self-consistent  time-dependent  analytical  theory,  it 
seems  to  us  that  the  computer  simulation  model  is  an  indispensable  tool  for  the 
explanation  and  guidance  of  future  CIA  experiments.  Our  model  can  also  be  used 
to  explain  the  fast  ion  energy  spectrum  observed  in  laser  produced  plasma  experi¬ 
ments.  The  fast  ions  are  a  result  of  collective  acceleration  by  the  high  energy 
electrons  emerging  from  the  laser  absorption  layer.  The  model  can  be  used  to 
relate  the  fast  ion  spectrum  to  the  electron  spectrum  for  diagnostic  purposes. 
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Fig.  1 


Fig.  2 


Fig.  3 


Fig.  A 


Fig.  5 


Fig.  6 


Fig.  7 
(a-e) 


FIGURE  CAPTIONS 


:  A  schematic  drawing  of  the  simulated  system;  $  is  the  electrostatic 
potential,  is  the  electrostatic  field,  d  is  the  length  of  the  drift 
tube;  the  plasma  source  can  be  a  proton  rich  thin  foil  or  a  puffed-in 
gas,  ionized  by  the  e-beam  or  by  laser  light. 

:  Flow  chart  of  the  computer  code  simulating  the  collective  ion  acceleration 
by  linear  electron  beams.  Space  charge  limited  emission  of  ions  is  assumed. 

:  The  configuration  of  the  electrostatic  potential  well  at  several  times, 
for  the  case  with  no  ions  present.  Initially,  there  is  a  transient  over¬ 
shooting  below  V  .  Then,  the  depth  of  the  potential  well  stabilizes  and 

2 

oscillates  around  the  beam  electron  energy,  Vq.  =  10  KA/cm  ,  Vq  ■  1  MeV. 

:  The  time  variation  of  the  virtual  cathode  amplitude  ( |  d>min ( t )  | )  for  the 
case  with  no  ions  present.  The  potential  stabilizes  after  approximately 
an  electron  crossing  time  along  the  drift  tube. 

:  The  time-dependent  behavior  of  the  various  current  components  in  the 

case  with  ions  present.  The  notations  are  explained  in  the  text.  The 

2 

injected  e-beam  current  is  10  KA/cm  . 

:  The  time  variation  of  the  ion  energy  distribution  function  in  four  energy 
ranges:  0-1,  1-2,  2-3  and  3-4  MeV.  The  maximum  energy  observed  is  A  MeV. 
The  gradual  transfer  of  energy  from  slower  to  faster  ions  is  marked  by 
sets  of  arrows.  The  units  are  logarithmic,  arbitrary,  but  relatively 
correct. 

:  The  configuration  of  the  electrostatic  potential  at  intervals  of  150  psec, 


for  the  case  with  ions  present.  A  conspicuous  feature  is  the  occasional 
appearance  of  short-lived  regions  of  positive  potential  ("virtual  anodes") , 
which  are  partially  responsible  for  the  ion  acceleration. 

Time  derivatives  of  the  ion  energy  spectra  in  several  energy  ranges. 

The  arrows  show  the  gradual  energy  transfer  from  one  range  to  the  next. 


FLOW  CHART  OF  THE  CODE  (in  C.G.S.  Units) 
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